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Abstract —The problem of base station cooperation has re¬ 
cently been set within the framework of Stochastic Geometry. 
Existing works consider that a user dynamically chooses the 
set of stations that cooperate for his/her service. However, this 
assumption often does not hold. Cooperation groups could be pre¬ 
defined and static, with nodes connected by fixed infrastructure. 
To analyse such a potential network, in this work we propose a 
grouping method based on proximity. It is a variation of the 
so called Nearest Neighbour Model. We restrict ourselves to 
the simplest case where only singles and pairs of base stations 
are allowed to be formed. For this, two new point processes 
are defined from the dependent thinning of a Poisson Point 
Process, one for the singles and one for the pairs. Structural 
characteristics for the two are provided, including their density, 
Voronoi surface, nearest neighbour, empty space and J-function. 
We further make use of these results to analyse their interference 
fields and give explicit formulas to their expected value and their 
Laplace transform. The results constitute a novel toolbox towards 
the performance evaluation of networks with static cooperation. 

Keywords — Cooperation; Static groups; Poisson cellular net¬ 
work; Thinning; Interference 

I. Introduction 

Cooperation between base stations (BSs) is receiving in 
recent years a lot of attention, due to its potential to improve 
coverage and spectral efficiency. It has shown considerable 
benefits especially for cell edge users that suffer from inter-cell 
interference. It is also expected to play a significant role due 
to the coming densification of networks with HetNets [1]. The 
concept of cooperation in the downlink implies that two or 
more BSs exchange user state information and data to offer a 
stronger beneficial signal with reduced interference. The total 
benefit depends on the amount of information exchanged, but 
also very importantly on the number and positions of nodes 
that take part in the cooperation. 

Recent studies have approached the problem of downlink 
cooperation with the theory of Point Processes [2] and Stochas¬ 
tic Geometry [3], where the network topologies follow a cer¬ 
tain probability distribution. Often, BS positions are modelled 
by a Poisson Point Process (PPP) with some fixed density over 
the entire plane. Specifically, Baccelli and Giovanidis [4] have 
studied cooperation between pairs of BSs. Evaluation for any 
number of cooperating stations has been done by Nigam et al 
in [5], Tanbourgi et al in [6] and Blaszczyszyn and Keeler in 
[7], In all these works, the common ground is the use of PPPs 
and the fact that the cooperation is driven by the user, who 
defines the set of stations for his/her service. 

However, the last assumption is not very realistic since 
it overburdens the backhaul/control channel with intensive 


communication between BSs. Furthermore, it is not very clear 
which user is served by which station. A simpler and more 
pragmatic approach is to define a-priori static groups of BSs 
(that do not change over time). In each group, the BSs may 
reliably communicate with each other by reservation of control 
channel bandwidth or installation of optical fibres between 
them. If the criterion for grouping relates to geographic prox¬ 
imity, BSs in a small distance will coordinate fast and will 
share a planar area of common interest. The important question 
is how exactly should these groups be defined? 

The idea is not new and suggestions have already appeared 
by Papadogiannis et al [8], Giovanidis et al [9], Akoum and 
Heath [10] and Pappas and Kountouris [11], Static clusters 
have also been considered for the uplink in [12], [13] and 
[14], In these works however, groups are formed neither 
systematically, nor optimally. Other works model the problem 
of dynamic clustering as a coalition game [15], 

We propose in our work a criterion for BS grouping that 
only depends on geometry. It is a variation of the Nearest 
Neighbour Model for point processes suggested by Haggstrom 
and Meester [16]. According to this, two BSs belong to the 
same group if one of the two is the nearest neighbour of the 
other. This assumption is reasonable for the telecommunication 
networks because it forms groups based on proximity. The 
variation we consider here limits the maximum number of 
elements K that can group together. Since the general problem 
is very complicated, after presenting the criterion in Section 
II, we focus on the case K = 2, where the BSs can be 
either single or cooperate in pair with another BS. This is 
already interesting, as it raises all the important questions of 
the general problem. It leads to fundamental results and is very 
challenging to study. 

In Section III we formally give the notions of single BSs 
and pairs and define the two point processes d/W and 
that result from a dependent thinning of the original process 
<1». The total interference of the network results from the sum 
of the interferences of the individual processes and <1>^ 2 K 
In Section IV we show that these processes are not PPPs 
and provide many of their structural properties: the average 
proportion of atoms from $ that belong to ‘I>'■ 1 and 4>' 21 , 
the average proportion of Voronoi surface related to each of 
them, as well as properties concerning repulsion/attraction. 
We continue in Section V with the interference analysis 
and obtain explicit expressions for the expected value of the 
interference created by each one of the two point processes. 
We also provide their Laplace Transform (LT) when they are 
constrained within a finite subset of A C K 2 . Finally, Section 
V concludes our work. All proofs of theorems can be found 


in the Appendix. 


II. Organising Base Stations into Groups 

Let us consider a Poisson Point Process (PPP) >I> in R 2 
with density A > 0. A realisation of the process can be 
described by the infinite set of (enumerated) atoms {xi}. 
Each realisation represents a possible deployment of single 
antenna Base Stations (BSs) on the plane. We wish to organise 
these BSs (or atoms) into cooperative groups C m (</>), with 
possibly different sizes, where size refers to group cardinality 
card (C m ). The index m enumerates the formed groups. We 
consider groups of atoms whose union exhausts the infinite set 
(f> and they are disjoint 
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We aim to find groups that are invariable in size and 
elements with respect to the random parameters of the telecom¬ 
munication network (e.g. fading, shadowing or user positions). 
In this sense, we look for a criterion that aims at network- 
defined, static clusters that differ from the user-driven selection 
of previous works. For this reason, we use rules that depend 
only on geometry. Based on these, an atom x £ f takes part in 
a group, based solely on its relative distance to the rest of the 
atoms <j> \ {x}. This geometric criterion is related to the path- 
loss factor of the channel gain. When a user lies at a planar 
point z and is served by BS x, the gain is equal to 

h (z,x) = a (z,x) d(z,x)~ /3 , z £ M 2 , x £ <j>, (3) 

where d(z, x) := \z — x\ is the Euclidean distance between the 
user and the BS, d (z, x) 13 is the path loss with exponent f3 > 2 
and v (z, x) models the power of channel fading. Both h and 
v hence refer to power, while the channel fading component 
is a complex number equal to (e.g. for Rayleigh fading 

v is exponentially distributed). 

A. The Nearest Neighbour Model. 

In our work we investigate grouping decisions based on a 
model proposed by Haggstrom and Meester [16], and further 
analysed in [17], [18], [19], the so called Nearest Neighbour 
Model (NNM). Given the realization we connect each atom x 
to its geometrically Nearest Neighbour by an undirected edge. 
This results in a graph Gnn, which is well defined, because 
for a PPP no two inter-atom distances are the same a.s. and 
hence each atom has a unique first neighbour. However, an 
atom can be the nearest neighbour for a set of atoms (possibly 
empty). 

The NNM has certain properties that make it a good candi¬ 
date for our purpose. (P. 1) The group formation is independent 
of the PPP density A. (P.2) The graph Gnn is disconnected, 
i.e. there always exist two atoms not connected by any path. 
(P.3) Each resulting cluster C does not contain cycles , it is a 
tree and hence the graph Gnn is a forest. (P.4) The graph 
contains a.s. no infinite component, i.e. it does not percolate 
[16, Th.2.1 and Th.5.2]. Consequently, the cardinality of each 
cluster is a.s. finite. (P.5) All atoms necessarily have a nearest 
neighbour. An example of a Gnn for a realisation of a PPP 
within a fixed window is shown in Fig.l. 
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Fig. 1: Example of cooperation groups using the Nearest- 
Neighbour Model 

Model Variation: Although the NNM guarantees that the 
groups are finite, there is no upper bound on the group size. 
When referring to a telecommunications network, however, it 
is often more reasonable to bound the maximum group size by 
a number K. For I\ = 2 this means that there can appear only 
single atoms and pairs, for I\ = 3 singles, pairs and triplets etc. 
It is this modification of the NNM that we propose and analyse 
in this work, because it will lead to more natural grouping of 
BSs. Algorithmically, starting by a realisation of single atoms, 
we first search for possible groups of two (pairs). These are 
formed when two atoms are mutually nearest neighbours of 
each other and we result in the case K = 2. From this last 
constellation, to shift to the case K = 3, we will consider the 
existing pairs. For each pair we will search among the single 
atoms to find if there exist some of them that have one of the 
pair as nearest neighbour. If so, we choose the geographically 
closest with this property and create a group of three. Observe 
that not all of the pairs will become triplets. We iterate in this 
way for larger I\. 

III. The special case of NNM groups with K = 2. 

From this point on, the paper will be devoted to the study 
of the simplest of static cooperation cases, the case I\ = 2. 

A. Singles and Pairs. 

The following definitions may apply to any process $. In 
our work we will analyse the case where $ is a PPP. For two 
different atoms x, y £ </>, if x is in Nearest Neighbour Relation 
(NNR) with y, that is, if 

y = argmin d(x, z), 

ze<f>\{x} 

we write x A y. If this is not true we write x y. We use $ 
instead of <j>, when we consider the whole set of realisations 
Cl (pj £ fl, <f> (w) = <j>), i.e. when we calculate probabilities. 
We will omit the dependence on $ or when it is clear from 
the context. 

Definition 1. Two atoms x,y £ <f>, are in Mutually Nearest 
Neighbor Relation (MNNR) if, and only if, x y and y % x. 



We denote this by x y. The two atoms then form a pair. 
In telecommunication terms we say that BSs x and y are in 
cooperation. 

Definition 2. An atom x is called single if it is not in MNNR 
(does not cooperate) with any other atom in </>. (Formally if 

for every y G 4>\ {a;} such that x \ y, it holds that y ■$* x). 
We denote this atom by x^. We then say that BS x transmits 
individually. 

In geometric terms, the NNR x A y holds if and only 
if (iff) there exists a disc B(x,r) with radius r = \x — y\ 
empty of atoms in 0. The nearest neighbour of x is unique 
in the case of PPP, because the probability of finding more 
than one atom on the circumference is zero. Furthermore, 
the MNNR holds, iff a symmetric relation is true i.e. the 
disc B (y, r) centred at y with the same radius is empty. 
We conclude that the MNNR holds for x,y, iff the area 
C ( x , y) := B ( x , |x — y\) U B (y , \x — y\) is empty of atoms. 
Its surface S is equal to S (C(x,y)) = n \x — y\ 2 (2 — 7). 
Here, 7 := § — ^ ~ 0.391 [18] is a constant number equal 
to the surface, divided by 7 r, of the intersection of two discs 
with unit radius and centres lying on the circumference of 
each other. An illustration of the above explanations is given 
in figure 2(a). On the other hand, figure 2(b) illustrates the 

single atom definition, where for any y for which x y, the 
disc with centre y and radius \x — y\ contains at least one atom 
of the process. 


(a) The atoms x and y are in MNNR (b) The atom in is a single and we 
x&y (i.e. x -X y, and y -X x). write w # ( Le - w-r y, but y ^ w). 

Fig. 2: Illustration of a pair of atoms and a single atom. 

With the above, and the empty space function for PPPs [3], 
we find the probability of two atoms being in pair. 

Lemma 1. Given a PPP $ with density A and conditioned 
that at x,y G R 2 there are tw’o of its atoms, the probability 
that these are in MNNR is equal to 

P (x&y) = t (4) 

where 7 := f - # 

Next, we give the global probability of any atom of $ 
being single or in a pair. This result is shown from the point of 
view of an atom at x, i.e. involving the Palm measure P x (see 
Appendix). We can heuristically understand the Palm measure 
as P x (-) = P(-|x G $) [20]. This description explains well 
its use, but the probability of {x G $} is always null. Still 


we can assume that there exists a point of 4> within a small 
neighbourhood around x (say a ball with radius e > 0 and 
centre x) and then calculate the required probability, taking 
the limit e —> 0. 

Theorem 1. Given a PPP $ with density A and x G $, there 
exists a constant p*, independent of A and x, such that 

P (x G> y, for some y G = P* ■ 

P(x|) = l -p*. 

Specifically, p* = ss 0.6215. 

We should remark that Theorem 1 gives also the percentage 
of points that are singles or in pair, for any planar area ACK 2 . 
This is independent of the point of view of an atom at x. It 
means that, in average, a percentage 37.85% of atoms of An$ 
are singles and 62.15% of atoms are in pair. This has been 
verified by Monte Carlo simulations over a finite (but large 
enough) window A. It can also be analytically evaluated, using 
/O) = i{ xeA} i n the expression (15), later on. The MNNR 
criterion leads, hence, to a reasonable splitting of the initial 
process into two processes, one with singles and another one 
with pairs of cooperating BSs. 

Now, we can define two new point processes <l> rl * and 
that result from the dependent thinning of the PPP $ with the 
MNNR criterion, using Definitions 1 and 2 

<b (1) = {x e $ & x is single}, 

<fA 2 ) = {x G $ & x cooperates with another element of $}. 
B. Voronoi Cells. 

It follows naturally to investigate the size of Voronoi cells 
associated with single atoms or pairs. A Voronoi cell of atom 
x G is defined to be the geometric locus of all planar points 
z G R 2 closer to this atom than to any other atom of ( f > 
[21], In a wireless network the Voronoi cell is important when 
answering the question, which users should be associated with 
which station. Let z rx <j>V) (resp. z rx fi^ 2 ' 1 ) denote the event 
that z belongs to the Voronoi cell of some atom of tfM'l (resp. 
0 <2> ). For the probability of these events we have analytical 
forms but not their numerical solutions, due to numerical issues 
related to integration over multiple overlapping circles. 

Numerical Result 1. The average surface proportion of 
Voronoi cells associated with single atoms and that associated 
with pairs of atoms is independent of the parameter A. By 
Monte Carlo simulations, we find these values equal to 


p 1 

(z^$ (1) ) 

« 0.4602. 

(6) 

p 1 

[z rx 

« 0.5398. 

(7) 


Interestingly, although the ratio of singles to pairs 
is 0.3785/0.6215 the ratio of associated surface is 
0.4602/0.5398, implying that the typical Voronoi cell 
of a single atom is larger that that of an atom from a pair. 
The last remark gives first intuition that there is attraction 
between the cooperating atoms and a repulsion among the 
single atoms. We will further analyse these observations in 
the next section. We close this paragraph giving an example 
of the association cells for the singles and the pairs in Fig. 3. 







Fig. 3: Example of planar association areas for singles and 
pairs for a realisation of a PPP in a square window of size 
t = 5. It is based on the Voronoi tessellation. 


For the pairs we consider the union of Voronoi cells of their 
individual atoms, as total association cell. 

IV. Characteristics of and <44 

In this section we consider each one of the two newly 
defined processes 4 1 ) and 4 2 ' separately and analyse their 
behaviour. Specifically, we try to establish possible similarities 
related to PPPs and for each one discuss issues of repulsion 
and attraction between their atoms. 


a typical atom of the process to the nearest other atom of the 
process [20]. We first analyse 4 2 ) and for this we denote by 
p( 2 ) ,x , and G^(r) the associated Palm probability measure 
and NN function of the distance r > 0, respectively. 

Theorem 2. The NN function of 4 2 ) is equal to 

G ( 2 ) (r) = p(2 )>*(d(x, $ (2) \ {z}) < r) 

= i_ e -^ 2 (2-7) ; ( 8) 

where 7 is the same constant as in Lemma 1. 

Hence, the NN random variable (r.v.) is Rayleigh dis¬ 
tributed, with scale parameter a = (2Xn(2 — 7 )). -1 ' 2 . Un¬ 
fortunately, we have not been able to provide a closed analytic 
expression for the NN function of the process <l> - 1 *, as we did 
for 4 2) in equation (8). In Figures 4(a) and 4(d) we plot the 
NN function of <4-0 using Monte Carlo simulations and '\ >!2> 
from the analytic expression in (8), respectively. Both plots are 
compared to the NN functions of the independent PPPs <4 11 
and <44 These are equal to (for i £ {1, 2}) 

G«(r) = P X (d(x, \ {a;}) < r) 

= 1 — e _Ai7rr ’ 2 . (9) 

To get an intuition how different the results are, observe that 
the exponent in the PPP case G t2> depends on the density 
A 2 = A/(2 — 7 ), whereas the expression for G 4 is the same 
but with exponent A (2 — 7 ). This observation explains why 
the use of independent thinning would be an inappropriate 
approximation in our case. 


A. Non-Poissonian behaviour. 

A first property is that both processes are homogeneous. 
This is due to the homogeneity of $ and because by definition, 
4 2 ) depends only on the distance between elements of >I>. 
Similarly for <44 Since the two processes result from a 
dependent thinning of a PPP $ they may not be PPPs. In 
fact we can show that they are not: Suppose that 4 2 ) is a 
homogeneous PPP. As shown in Theorem 1, the percentage 
of its atoms in MNNR with some other atom of 4 2 ) should 
be around 62.15%. However, by definition, 100% of the 
elements of ‘l>' 2> are in MNNR and we conclude that < 4 2 ) 
is not a PPP. For 4 1 ) the argumentation is not as simple. 
Nevertheless, we can show this with Monte Carlo simulations 
on the average number of single points, which is far from 
the number 1 — p* = 37.85%, and also using the Kolmogorov- 
Smirnov test [22], which shows that the number of 40 atoms 
within a finite window is not Poisson distributed. 

In what follows, we will use the PPP as reference process 
for comparison with the two new processes. For this, we 
consider two independent PPPs, 4 11 and 44 that result from 
independent thinning of the original PPP with probability 1— p* 
and p* respectively. This is motivated from Theorem 1. These 
PPPs will have density Ai = A(1 — p*) and A 2 = A p*, which 
gives the same average number of points as the processes 40 

and 44 

B. Nearest Neighbour function. 

The Nearest Neighbour function (NN), denoted by G, is 
the cumulative distribution function (cdf) of the distance from 


C. Empty Space function. 

The empty space function (ES), denoted by F, is the cdf of 
the distance from a fixed planar point z £ R 2 to the nearest 
atom of the point process considered [20]. By stationarity of 
40 and 44 the function F does not depend on z, so we 
can consider the typical user (point) at the Cartesian origin. 
Unfortunately, we have not been able to derive analytical 
formulas for either of the two F lf 1 and 44 The ES functions 
of the independent PPPs are denoted by 4 1 j and F' 2 ' and they 
are equal to the expression of the NN functions in (9). This is 
a property of the PPPs [20], 

We have used Monte Carlo simulations to plot the ES of the 
two processes. In Fig. 4(b) we show the comparison between 
j4) and 40 , which - unexpectedly - seem very close to each 
other. The same comparison in Fig.4(e) shows also closeness 
of fit, although less tight, between F (2> and 44 Since we 
have no analytic expressions, we are tempted to consider the 
ES function of the independently thinned PPPs as a reasonable 
approximation for those of 4 1 > and 44 


D. The J function. 


The two functions NN and ES can be combined into a 
single expression known as J function. The latter is a tool 
introduced by van Fieshout and Baddeley [20] to measure 
repulsion and/or attraction between the atoms of a point 
process. It is defined as 


J(r) 


1 - G(r) 

1 - F(r) ’ 


GO) 














Fig. 4: NN function, ES function and J function for the processes (bW, and their comparison with the PPPs <l>• 1 ) and 2 h 


In the case of the uniform PPP, G (r) = F (r) and J (r) = 1, 
as a consequence of the fact that the reduced Campbell 
measure is identical to the original measure. Hence the J 
function quantifies the differences of any process with the PPP. 
When J(r) > 1, this is an indicator of repulsion between 
atoms, whereas J(r) < 1 indicates attraction. We plot in 
Fig.4(c) the / function of d?) 1 ) and in Fig. 4(f) that of $( 2 ). 
From the figures we conclude that exhibits repulsion for 
every r > 0, and d>( 2 ) attraction everywhere. However, note 
that the attraction in the case , 1> : 2 ' is due to the pairs formed. 
If we consider a new process having as elements the pairs of 
d> !2 / this process of pairs exhibits repulsion everywhere. 


the Cartesian origin due to stationarity. To keep our setting as 
general as possible we describe by use of measurable functions 
/ : R 2 — y R + , g : R 2 x R 2 —y R + the signals transmitted 
by a single BS or by a pair and received at the typical user. 
The sum over the entire (of each) process gives the random 
variables of interest 


j(t) 

X(2) 


/(*)» 


1 

2 


H 9{x,y)t 

x£0( 2 ) y£0( 2 ) 



(11) 

( 12 ) 


V. Interference Analysis 

The purpose of the previous analysis was to develop the 
tools necessary for use in a communications context. Within 
this context, the cooperating BSs will have a different influence 
on the interference seen by a user in the network, than those 
operating individually. The current section will focus on the 
interference field generated by d>W and d>( 2 ). As shown in 
Section IV, the two processes have a different behaviour 
compared to a PPP and this is why approximations based on 
independent thinning of $ will not bring accurate results. We 
thus have to resort to a more direct approach. 

We denote by X' 1 j and I { ' l ! . the interference field generated 
by <l> n ^ and respectively. The typical user is chosen at 


The 1/2 in front of the summation in (12) prevents us from 
considering a pair twice, whereas f implies that y £ <// 2 )\{:e}. 

We now give some practical examples for the functions 
under study. First of all, the function for the individual BS in 
this analysis will be equal to 

f(x) = h(0, x) (13) 

where h(z,x) was defined in (3). This is the signal from a 
single antenna BS that lies at distance |0 — a) from the user 0. 
The fading power is(0,x) follows the exp(l/P) distribution, 
where P is the BS transmission power [23], The cooperation 



























signal is more interesting. We can consider the following cases 


9{x,y) = 


h(0,x) + h{ 0,y) 
max{/i( 0 , x), h( 0 , y)} 
l on h(0,x) + (1 - lon)h(0,y) 

y/h( 0 , x)e l9:r + y/h( 0 , y)e l9y 


[NC] 

[OF1] 

[OF2] . (14) 
[PH] 


[NC] refers to the no cooperation case, where both BSs of the 
pair behave individually. [OF1] refers to the case where the BS 
with the strongest interfering signal is actively serving its user, 
while the other is off (alternatively we could replace max by 
min to consider the weakest signal of the two). [OF2] is again 
a scenario with one of the two BSs active and the other off, 
where the choice is made randomly by a r.v. with E [l on ] = q 
(e.g. q = 0.5 for fairness). Finally, [PH] is the case where the 
two complex signals are combined in phase as well [4]. Here, 
i is the complex unit and 0 X , 6 y are the uniformly random 
phases of the signals. Observe that the above signals can be 
generalised to include MIMO transmission as well. 


plane. Then by direct application of the expected value of a 
function y, we can write the LT as an infinite sum of terms 
y(n)P (./V(A) = n), where y(ri) is the value of the function 
when n atoms appear in A and P (N(A) = n) is the probability 
that this event occurs. 

The method can be seen as an approximation of the LT 
of and <fA 2 ). It is an open question to prove that when 
A —> R 2 the LT of the finite process converges to the LT we 
are looking for. We write d>M ~ <|>''', for A regular and large 
enough. However, since we are often interested in simulating 
and analysing only finite areas and also since the conjecture in 
the limiting case sounds reasonable, the result we present here 
has a great importance. It fully characterises the distribution 
of the interference (and is provable for finite windows). 

For a finite number n of known planar points Xi,... ,x n £ 
A (that are potentially occupied by atoms, when the latter are 
chosen uniformly within the area), we define the functions 

H M : (R 2 ) n -» (K+) n and J : (K 2 )” (K +) n2 


A. Expected value of 2A 1 ) and I'- 2 '. 

The next Theorem gives an exact integral expression to 
the expected value of the interference field generated by the 
singles and the pairs. The proof uses the Campbell-Little- 
Mecke formula. Lemma 1, and Theorem 1. 


Theorem 3. The expected value of the interference field 
generated by T4 1 ) anc j q>( 2 ) i s given by 


jW 

j(2) 


= (1 ~P*) / E[f(x)]Xdx, 


(15) 


/ R 2 


E [g{x, y)} e-^-^^-^XdyXdx. 

(16) 


The expected value can be finite or infinite, depending on 
the choice of f(x) and g(x,y). Observe that for [NC] and 
[PH] the expected interference has the same value. 


H\ n \xi, ...,x n ) = 


if Xi satisfies the ’’single” relation, 
otherwise. 


The ’’single” relation is satisfied in the deterministic case, ex¬ 
actly as in Definition 2 if, for every Xj 7 ^ Xi for which ay Xj 


it holds that Xj Xi. Then H^ = (H 


(n) 


,hY) t . 


Ji n \x 1 , 


J( 0 ... 1 ...), if Xi is in ’’pair” with Xj. 
\( 0 ... 0 ), in other case. 


The ’’pair” relation is satisfied in the deterministic case, as 
in Definition 1 if a;, f> Xj , for some i 7 ^ j. Then = 
(j[ n \ , jY) T , where J T is the transpose vector of J. 


Theorem 4 (Laplace transform). Consider a PPP $ with 
density X, a subset A C M 2 and the functions f, g (e.g. (13) and 
(14)) related to the single atoms and the pairs respectively. Let 
F^(x = {f{xi),...,f(x n )), G[ n \xi,... ,x n ) = 

(g(xi,x 1 )... g(xi, x n )), and G (n) = (G^ n) ... G^). 


Corollary 1 (Intensity measure). The intensity measures for 
and denoted by and M^ respectively, are 

equal to 

M^\dx) = (1 — p*)Xdx, 

M^ 2 \dx) = p*Xdx. 

B. Laplace transform of and 1^. 


(18) 


The LT of the interference for is equal to 


E 

—sX^ 

e s± 

= e ~ xs(A) f 1 + A / E 

e sf(x) 

(17) 


\ J A 

■ 


, A 2 
dx + Y - 


00 \n r 


E 


sF^.H^ 


dx 1 ... dx r 


As a final result we present our findings related to the 
Laplace transform (LT) of the interference from <I > (1 '> and 
We derive here the exact LT of the interference when 
considering a finite subset A C M 2 (window) and the related 
random point measure 

4>y 1 = {single atoms of $(A) = 4> n A inside A.} 

= {atoms of $(A) in MNNR with another atom in A.} 

A sketch of proof is given due to space constraints. For finite 
subsets, the atoms of the PPP 4> are distributed i.i.d. uniformly 
within A. We condition on the number of atoms that appear 
within A, and make use of the fact that is a PPP (or some 
other process with known counting measure). We also consider 
the MNNR only among the atoms in A and not on the entire 


The LT of the interference T^ 2 ' 1 for is equal to 
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We finish this section by evaluating the expected value of 
the interference using Theorem 3 and the proposed expressions 
for f(x) in (13) and for g(x, y) in (14). We further compare 
the results from the numerical integration with those by Monte 























Nevertheless, the two scenarios do not numerically defer much 
for /3 = 4, as shown in the figure. 

VI. General case and Conclusions 

The analysis for the dependent thinning of a PPP using the 
NNM can be generalised to include groups of atoms with size 
greater than two. For the computation of such probabilities, 
certain difficulties are raised due to the overlapping of more 
than two discs with different radii. However, the methodology 
remains the same. Results on the percentage of atoms, Voronoi 
surface and repulsion/attraction, which are independent of the 
PPP density, can be derived by Monte Carlo simulations. The 
same analysis for the expectation and the LT of the interference 
can also be followed in cases of larger groups. 


Fig. 5: Interference generated by the single atoms outside a 
ball of radius R. The upper plot shows the case with path-loss 
exponent 3 = 2.5 and the lower plot the case with 3 = 4. 



Fig. 6: Interference generated by pairs of atoms outside a ball 
of radius R. The upper two plot shows the case with path-loss 
exponent /3 = 2.5 and the lower two plots the case with 3 = 4. 
We show the two cooperation scenarios [NC] and [OF1], 


Carlo simulations within a finite - but large enough - window. 
The chosen density is A = 0.1 [atoms/m 2 ] and the window is 
a square of size 100 x 100 [m 2 ]. 

Specifically for I \, we first transform in polar coordinates 
(15) and then perform numerical integration, for the function 
f(r x ) = with i/(0, (r x , 0 X )) := 1. Here, R 

is a positive distance, within the interval R £ [0.5,5]. The 
indicator function is used in order to calculate the interference 
created by singles outside a ball centred at 0 and of radius 
R. The evaluation is shown in Fig. 5 with a continuous line. 
The results from the simulations are shown by the star-dotted 
curve. With the appropriate window, the expression in (3) gives 
almost identical results with the simulations. 

Similarly for 1-2, numerical evaluation and simulation re¬ 
sults using g{r x ,r y )l[ r;ctry> R} are shown in Fig. 6. For the 
numerical integration, we evaluate the two cases [NC] and 
[OF1] for different path-loss exponent. Obviously interference 
from [OF1] is always less than [NC] since it is received only 
from one of the two BSs of each pair, while the other is silent. 


Altogether, in this paper we proposed a method to define 
static groups of cooperating Base Stations based on the Nearest 
Neighbour Model and analyse them with the use of stochastic 
geometry. Many structural characteristics were derived for the 
case of singles and pairs. These can be further used in the 
analysis of SINR models and their performance evaluation. 
We gave some first results by calculating the expected value 
and the Laplace transform of the interference from different 
groups. Further analysis and applications in this direction is 
the current research work of the authors. 
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Appendix 


Proof of Theorem 1. 

Let us condition on the fact that an atom lies at the planar 
point x and then find the probability that this is in pair (single). 
This is written as P x ^a: £> y , for some y £ 4> \ {a;}^, where 
P ;,: is the Palm measure. We use Slivnyak-Mecke’s theorem [3] 
, which states that <1> \ {x} under P x has the same distribution 
as $ under P. 


for some y £ $ ) =E (^xAy, for some ye *} 
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0.6215. 



















Campbell-Little-Mecke formula, and Lemma 1, to find that 


Equality (a) holds because for PPPs the nearest neighbour of an 
atom is a.s. unique, so we cannot have three atoms x, y, w £ <j> 

with y f w, such that x A y and x A w. Equality (b) comes 
from Campbell’s formula [3] and in (c) we use Lemma 1. 

Similarly, to find the probability of x £ $ being single 

P 35 =1 — P x (x -O- y, for some y \ {a:} € 

1 — 7 

= --- « 0.3785. 
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Palm measures and the NN function of A 2 ). 

Let (fljJ 7 , P) be the underlying probability space and $ 
a PPP The Palm measure of $ is P x . We denote by pW’* 
and p( 2 )’ x the Palm measure of A 1 * and A 2 ), respectively. 
Let us first consider P^ 2 A X and use the heuristic definition 
p( 2 )’ x '(.) = P(.|a; £ A 2 )), where we condition on the existence 
of one atom of $ in a small neighbourhood of x. By Definition 
1, x £ </>( 2 ) if, and only if, A holds, where 


A x = {there exists some y £ 4>\{x} such that x ■£> y} 


E 

ye<l>\{x} 


1 , 4> , 

{x-H-y} 



Also, from Theorem 1, P X (A) = P* for PPPs. For every 

r 6 j, 

p( 2 )> x (r) = p( r \ x g $,A) 

= P(T, AN g $) 

P(A]ze$) 

= p^r, A) (19) 

P X (A) 

= p x (r|A) 

= p x (r, a) 4 - 

For a rigorous proof of the Palm measure for <!>' 2) , we refer 
the reader to the last part of the Appendix. 

To calculate the NN function of A 2 ), let us take r > 0 and 
consider T = {d(x, A 2 )\{x}) < r}. Then, 
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where the last equality is due to the fact that p* = • 


We further remind the reader that x £ (f 1; if, and only if, 
the event B x holds, where 


B x = { for every y £ </>\{x} such that x —> y, y x}. 

It is possible to give a similar expression to pW’ x in the same 
way as we did for p( 2 )’ x . For every T £ T, 


p( 1 )- x (r) =p x (t\b x ) 


P X (T,B X ) 
(1 ~P*) ' 


It has not been possible however to make a similar analysis to 
get the NN function of A 1 ). The reason is that it is not easy 
to precisely define the nearest neighbour atom in this case. A 
property similar to that in equation ( 20 ), which was crucial in 
the previous proof, could not be found. 


G ( 2 ) (r) = P( 2 )’ x (d(x,$ ( 2 ) \{a;}) < r) 

'= P x (d(x,^\{x})<r,A x )^ 

= W E 1 {^,‘& (2) \W)<^^!/} ) pp- 

\y£<5>\{x} ) 

For every <f> and y £ (j)\{x}, if x ■£> y, it is true that 

d{x, </> ( 2 ) \{ 4 ) = d{x, AM) = d(x, y). ( 20 ) 


Proof of Theorem 3. 

Let us start with T^\ We observe that 

I(1) = E /(4 = E^) 1 ^} p °- s -> 

xgi i>W 

and that, for every x £ <f> 
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We use this observation, Slivnyak-Mecke’s theorem, the 







By the reduced Campbell-Little-Mecke formula and Slivnyak- 
Mecke’s Theorem, 
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J R 2 

= (1-P*) [ E[f(x)]Xdx, 
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where (a) and (b) come from the proof of Theorem 1 . 
For JW, we make the observation that 

* * 
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Then, as previously (and using Theorem 1) we can calculate 
the expected value of the interference from pairs 
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(c) uses Campbell’s formula and (d) comes from Lemma 1. 


Proof of Corollary 1. 

Let us take A G £>(R 2 ). We use in (11) the function 
f(x) = which indicates whether the atom x belongs 

to A or not. We thus have 


and we conclude that 

M^(dx) = (1 — p*)Xdx. 

For we consider again the function g{x , y) = l{ x gA}- 
Given that, for every x G 0 (2) , £, 6 </ ,( 2 )\{x} 1 (2) = L 

{x y} 

equation ( 12 ) takes the form 

I (2) = 1{*GA} E 1 r AV 

:rE<j!>( 2 ) y£(j)( 2 )\{x} ^ 

= E ^A}- 

a?E<j!>( 2 ) 

Its expected value is flip’s intensity measure 
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The righthand side in (16) of Theorem 3 is 


f f g(x,y)e- x ^ x - y \ 2 ( 2 ~^XdyXdx = 
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[ [ t^e-^-^ 2 ^-^XdyXdx = 

J A J R 2 

[ [ e-^-yf^ 2 -^XdyXdx = [ p*Xdx. 

J A J R 2 J A 


We conclude that 


M^ 2 \dx) = p*Xdx. 

This result provides an analytical argument of the stationarity 
of the processes <I» ( 11 and <I> f 21 through an explicit formula of 
their intensity measure and their intensity coefficient. More¬ 
over, using the intensity measure, we can find the average 
number of points of <t> ■ 1 ) or over any Borel set. This 
clarifies our discussion at the end of section III.A under 
Theorem 1. 


2 (1) = £ f{x) 

ccE^ 1 ) 

= E 

which counts the number of elements of 01 1 ' inside A , and, 
hence, its expected value is $W’s intensity measure 
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E 1 1^6A} 
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$ (1) (A) 


M {1 \A). 


The righthand side in (15) of Theorem 3 is 



p*)Xdx 


[ (1 — p*)Xdx , 

J A 


Palm Measure for - rigorous proof 

In this section, $ is a stationary point process, with inten¬ 
sity A > 0, and Palm measure given by P°. For a realisation 
f = {x n ) and a fixed x G R 2 , we denote by (j> x := {x n + x}, 
the translation of <j> by x. 

In the previous sections, we have defined the processes 
<I» ! 1 ' and $ 11; only by geometric means and not making use 
of the probability law that governs <l». What we have said is 
that given a realisation <f>, 

cj)'- 1 ^ = {x G f | x is single }, 

0 ( 2 > = {x G f | x cooperates with another element of <fi}. 

Denoting by M the space of realisations, it is possible to 
consider 


0 HA 0 ( 1 ) 0 HA 0 < 2 > 
























as measurable mappings M —► M. Since we use mappings 
over the space of realisations, we can not say much about the 
probability laws governing ( 1 > f 11 and < 1 > ! 2 \ which we denote by 
pt 1 ' 1 and p( 2 \ respectively. Let us consider the event Tg J. 
Each pW is a measure 

P (i) :T ->■ [0,1] 

J (21) 

r -> p (i) (r) = P($ w e r). 

We define rW g Jas := {<p \ </>(*) £ T}, therefore, 

P W(D = P G r«) = E (l { * e rw } ) ■ (22) 

(Let us note here that this last equation gives us a simple way 
to approximate pM using the Monte Carlo method.) 

Since $ is stationary, the Palm probability of $ is defined 
as follows [24, p.l 19], [20, p.51,(24)]. For every A £ B(] R 2 ), 
and every fg J, 

E Y l{x 6 A}l{*_. 6 f} = A5(A)P° (f) . (23) 

Using the above definition, we want to find a probability 
measure p( 2 )>° : p —>• [0,1] such that, for every A £ B( R 2 ), 
and every V £ T, 

E Y M^A}t { ^ ler} =p*XS(A)P^ 0 (T). (24) 

(Observe that we use the different events T and E for 
reasons that will be clear later in the proof). Let us denote by 

Aq := {<f> | there exists y £ 4> such that 0 A y}. For a given 
realisation <fi, it holds x £ q^ if, and only if, </>_ x £ Aq. 

We assume that the Palm measure for is p( 2 )’°(r) := 
P° (r( 2 )|„4 0 ) and we want to verify that the measure satisfies 
(24). Since, P° (_4o) = P* >0, it is definitely a probability 
measure. Let us take A £ £>(R 2 ) and r g T. We start by the 
righthand side of (24) 

p*A5(^)P (2) ’°(r) = p* XS (.A) — 

= XS(A)P° (t( 2 \a) 

= A5(A)P° (f) 

= e^ 

= E Y, V*e A }^{$_«er( 2 Mo} 

= E Y 

:r6<E> 

= E Y, 1{ X <=A}1 j^^grt 2 )} 

- E H l{xeA}l| $ (2) 6r |. 
xe $( 2 > 1 * J 

So we reached the lefthand side of (24). In the above, (a) comes 
by replacing F := {r( 2 \Alo}, (b) because x £ (/A 2 ) if, and 
only if (j>_ x £ Aq, and (c) because by definition </>_ x £ T^ 2 -* 
if, and only if <£ T. 


In a similar way, for p( 1 )>° : T -a [0,1], and for every T £ 
T we can show, using the same arguments, that p( 1 A°(r) := 
P° (r^ 1 )) Bq) is a Palm measure for f&W. 
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